clc;clear
x=randn(1,1000);
hist(x)
[mu,sigma] = normfit(x)
d=pdf('norm',x,mu,sigma);
figure
plot(x,d,'.')
